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Abstract. We study the entanglement evolution between two harmonic oscillators 
having different free frequencies each leaking into an independent bath. We use 
an exact solution valid in the weak coupling limit and in the short time non- 
Markovian regime. The reservoirs are identical and characterized by an Ohmic spectral 
distribution with Lorcnts-Drudc cut-off. This work is an extension of the case reported 
in [1] where the oscillators have the same free frequency. 
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1. Introduction 

Entanglement is considered a fundamental feature and a resource in the context 
of quantum computation and information physics [2J. However entanglement 
between quantum systems is easily destroyed when their unavoidable interaction 
with the external environment is taken into account [3]. Features and timing of 
the disentanglement phenomenon depend strongly on the particular physical system 
and environment under investigation. However, even if disentanglement cannot be 
avoided, quantum information and computation tasks can still be performed when 
the disentanglement time is much longer than the time needed to run the quantum 
task. From this point of view it becomes of fundamental importance to investigate 
mechanisms and features of disentanglement in those systems used in the context of 
quantum computation and information. 

Continuous variable (CV) quantum systems are some of the possible candidates 
for quantum protocols pE]. We consider here the problem of two independent harmonic 
oscillators each one leaking in a bosonic structured reservoir [5j El [7J El |9j [TUl [11]. This 
situation has been considered also by the present author in [lj. However here, as an 
extension, we study the case of oscillators with different frequencies, initially in an 
entangled TWB state. Our aim is to investigate the entanglement evolution in the case 
of two identical high-T Ohmic reservoir with a Lorentz-Drude cutoff. In this kind of 
systems, whose dynamics is determined by a generalized Hu-Paz-Zhang master equation 
[12] . all the initial Gaussian states remain Gaussian during the evolution. For this reason 
we are entitled to use the separability criteria for CV systems introduced by Simon [13J, 
as a marker of entanglement between the oscillators. 

The paper is organized as follows. In Section II we introduce the physical model 
and its master equation in the weak coupling limit. In Section III we review the concept 
of two-mode Gaussian squeezed state, or Twin Beam state, and provide the solution to 
the master equation. In Section IV we introduce the separability criteria for bipartite 
continuous variable systems and the Lorentz-Drude spectral function. Moreover we 
study the entanglement dynamics as a function of the initial amount of entanglement 
and the relative values of the free oscillator frequencies and the spectral function cut-off 
frequency. Finally in Sec. V we give a brief summary of our work. 

2. The master equation 

Our system consists of two non-interacting quantum harmonic oscillators with 
frequencies u\ and u 2 - Each oscillator is coupled to its own bosonic reservoir. The 
total Hamiltonian can be written as 

H = ^ fajOjaj + ^2 ^i^bjk + ^2 7i*( a i + a ])( b jk + b]k)> (!) 
j j,k j,k 

where j = 1,2, uji k and u 2k are the frequencies of the reservoirs modes, a, (a,-) and 
bjk (bj k ) are the annihilation (creation) operators of the system and reservoirs harmonic 
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oscillators, respectively, and 7^ describe the coupling between the j-th oscillator and the 
k-th mode of its environment. We assume reservoirs with the same spectral structure 
and equally coupled to the oscillators. The dynamics of the two oscillators can be 
described through the following non-Markovian local in time master equation [12] 



*-rMXliP(t)}-ijMXjAPi,P(t)}}, (2) 

where p(t) is the reduced density matrix, H® is the free Hamiltonian of the j-th oscillator, 
and Xj = ^(a>j + a ]), Pj = "^( a ] ~ a j) are t ne quadrature operators. The interaction 
with the reservoirs is taken into account through the time- dependent coefficients of Eq. 
fl2). The quantities Aj-(t) and Uj(t) describe diffusion processes, jj(t) is a damping term 
and Tj(t) renormalizes the free oscillator frequencies Uj. For environments in thermal 
equilibrium and in the weak coupling limit in (rj(t) being negligible), they read 

ft p+00 

Aj(t) = a 2 ds dujJ(u)[2N(u) + 1] cos(a;s) cos(u j-s), 
Jo Jo 

ft r+00 

Hj(t) = a 2 ds duoJ((jj)\2N(uj) + 1] cos(a;s) sm(ujs), 
Jo Jo 
ft r+00 

7j(t) = a 2 ds duJ(u)sin(us)sm(ujjs). (3) 
Jo Jo 



3. Gaussian states and separability condition 

As initial states for our system we consider the class of two-mode squeezed states (or 
twin-beams TWB) obtained applying a two-mode squeezing operator to the vacuum 
state of the oscillators [4J. They belong to the class of Gaussian states characterized by 
a Gaussian characteristic function Xo(A) = exp{ — |A T cr A}, where cr is the covariance 
matrix and A T = (xi,pi, x 2 ,P2)- For a TWB state we have 

with Ao = Diag(a,a), Co = Diag(c, — c), a = cosh(2r)/2 and c = sinh(2r)/2. The 
TWB state is thus determined only by the squeezing parameter r which also determines 
the initial amount of entanglement. The bigger is r, the larger is the initial entanglement. 
Using the characteristic function solution [TJ] , it is trivial to verify that the characteristic 
function at time t maintains its Gaussian character Xt{A) = exp{ — |A T cr t A}, where 



(i) 




with 



a t = ( C V 7m I > ( 5 ) 



At " Aoe H { -A« + n« A« + AS-ni? {<>l 
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( ce-( ri+r2 )cos[(u;i + u 2 )t\ ce~^ +r ^ sin[(wi + u 2 )t] . 
* ~ ^ - ce -( r i+ r 2) sinf^! + u 2 )t] ce~^ +r ^ cos[( Wl + u 2 )t\ 1 ' ' ' ' 

with rj(t) = 2 J Q *7(s)ds. Moreover, because we are interested in the short time non- 
Markovian dynamics only, we defined Ap (t) ~ J * Aj(s)c?s, and the following secular 
coefficients A^(i) — / * Ai(s) cos^ti^t — s)]ds, A*V(t) ~ J* Aj(s) sin[2o;j(t — s)]ds, 
11$ (t) ~ /* ni(s)cos[2a; i (t - s)]ds and IlS(t) ~ jj IL;(s) sin[2a;,(t - s)]ds. Details 
of the calculations can be found in [TJ [TTJ [13] . If a suitable environment spectrum is 
provided, all previous coefficients can be evaluated analytically in the high temperature 
limit. This is the case of the Lorentz-Drude Ohmic distribution we introduce in the 
next section. In equations (6) and (7) we omitted the explicit time dependence for all 
the appearing coefficients. 



4. Entanglement dynamics 

In this section we investigate the entanglement dynamics between the oscillators using 
the separability criterion of Simon [13]. This criterion is well-suited in the context of 
two-mode Gaussian states because it represents a necessary and sufficient condition for 
separability and it depends only on the analytic form of the covariance matrix. In the 
case of the time-dependent and non symmetric covariance matrix (jSJ), the Simon criteria 
is equivalent to the following algebraical inequality [5] 

4 (1) Aj 2) ] + '-- i^<- i^ 2 - tjaWtp.ta( 2 )' 



S{t) = det[A t w Af] + (- - | detC 4 |) 2 - Tr[Af>JCtJAl'JCjj] 



det[A 1; ] +det[A 

. > (8) 



with 




(9) 



S(t) is called separability function. When S(t) < the state is entangled, otherwise it 
is separable. 

We now consider a particular situation with two reservoirs in thermal equilibrium 
at high temperature T (hu c << ksT) characterized by a Lorentz-Drude Ohmic spectral 
function [3] 

J M = ^-^T772 ( 10 ) 



7T UJ Z + UJ 



c 



where u c is the cut-off frequency of the distribution. Hitherto we fix the temperature 
through the condition ksT/hw c = 100 and the coupling constant to the value a = 0.1. 
Moreover we use a dimensionless time variable r = u c t for the evolution of the 
separabilty. With these choices we have only two free parameters Xj = u c /ui (i = 1,2), 
namely resonance parameters, describing the relative positions of the free oscillator 
frequencies Ui with respect to the spectral cut-off frequency u c . In P has been observed 




Figure 1. (Colors online) Separability function 5 vs scaled time r = uj c t with r = 1 
and xi = 1. The three curves are relative to different values of x 2 . The blue solid 
line correspond to xi = 1, the red dashed line to £2 = 10 and the black dotted line to 
x 2 = 100. 



the existence of two different resonance parameters regions characterized by different 
qualitative and quantitative behavior of entanglement, x\ = x 2 > 1 and X\ = x 2 << 1. 
Exploiting this result we investigate three different regimes: 1) x±, X2 > 1, 2) x±, x 2 << 1 
and 3) x\ « l,x 2 > 1. 

A second degree of freedom for our analysis is the choice of the amount of 
entanglement in the initial TWB state, given by the value of squeezing parameter r. In 
principle we do not have any limitation in the choice of the value of r, however in real 
situations it is not possible to realize two- mode squeezed states with r > 2 [3]. On the 
other hand when r ~ 0.01 the amount of entanglement is so small that, at least in the 
high temperature limit, disentanglement is a very fast process almost independent from 
x\ and x 2 . For these reasons we restrict our investigation to the intermediate region of 
0.01 < r < 1. 

We start the analysis considering the case xi,x 2 > 1. We show the entanglement 
evolution in Fig. [1] where we fixed the squeezing parameter to the value r — 1. In this 
regime the dynamics is characterized by entanglement sudden death (ESD) appearing 
already in the short time non-Markovian region for any value of r. We also fixed the value 
of the parameter X\ = 1 and plotted three curves correspondent to x 2 = 1 (solid blue), 
X2 = 10 (dashed red) and x 2 = 100 (dotted black). We observe that the disentanglement 
time r dis , defined as Sfais) = 0, only slightly increases as the adimensional frequency 
detuning Ax = \xi — x 2 \ increases. 

The case xi,x 2 « 1 is reported in Fig. [2j As emphasized in detail in [T], this 
regime is characterized by a more variegate entanglement dynamics, showing oscillations, 
ESD and revivals. We report the situation relative to an initially small amount of 
entanglement correspondent to r = 0.04 and X\ = 0.1. The three curves correspond 
to x 2 = 0.1 (solid blue), x 2 = 0.2 (dashed red) and x 2 = 0.3 (dotted black). In 
all cases there are entanglement oscillations and ESD, while an entanglement revival 
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Figure 2. (Colors online) Separability function S vs scaled time r = u> c t with r = 0.04 
and X\ — 0.1. The three curves are relative to different values of x^ ■ The blue solid 
line correspond to X2 — 0.1, the red dashed line to X2 = 0.2 and the black dotted line 
to x 2 = 0.3. 




Figure 3. (Colors online) Separability function S vs scaled time r = uj c t with 
r = 0.1. Blue solid line correspond to the choice xi = x 2 = 0.1, red dashed line 
to xi = 0.1, X2 = 1, and black dotted line to x\ = x 2 = 1 



is also present when x 2 = 0.2. As in the previous regime of Fig. [T] we notice that 
disentanglement time is larger when Ax increases with fixed value of X\. However here 
the differences in the disentanglement time are much more evident. Also oscillations are 
damped out as the detuning increases. 

Finally we look at the regime x\ « 1 and X2 > 1 whose results are shown in 
Fig. |3]with the choice r = 0.1. The blue solid line corresponds to x\ = X2 = 0.1, the 
dotted black line describes the X\ — x-i — 1 case, and the red dashed line represents the 
intermediate case of x\ = 0.1 and x-i = 1- Here we can observe two completely different 
behaviors of S(t) relative to the cases x± = x 2 = 0.1 and X\ = x 2 = 1. The former case 
is characterized by evident entanglement oscillations and a much longer survival time 
of entanglement. Instead the latter case shows a fast ESD without signs of oscillations. 
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The intermediate case of x\ ^ x 2 possesses intermediate features: weak oscillations and 
a fast ESD. 

The three parameters regimes analyzed display all the possible features in the 
entanglement dynamics for our particular system: entanglement oscillations, ESD and 
revivals. Thus when the oscillators frequencies are detuned no new qualitative behavior 
emerges, compared to the equal frequency case [TJ. On the other hand we stress the 
existence of evident and non-negligible differences. This does not happen in the regime 
x\,X2 > 1, where the value of disentanglement time increases only slightly as the 
oscillator detuning Ax increases. 

Quantitative differences characterize the opposite case, i.e., the regime x±,x 2 « 
1. Even a slight detuning (Ax = 0.1) can change drastically the value of the 
disentanglement time and damp entanglement oscillations. Again this feature is present 
in Fig. |3] where these discrepancies are even more pronounced. 

From our investigation we can extract two main results. When the oscillators have 
different frequencies (ui ^ u 2 ), the entanglement dynamics shows a behavior which is 
intermediate between the cases of frequencies both equal to uj\ and both equal to oj 2 . 
Moreover the leading role in the disentanglement process seems to be played by the 
oscillator with higher resonance parameter x, which forces the entanglement to evolve 
in a way similar to the case of both parameters equal to the larger between the two, as 
shown in Fig. El 

5. Summary 

In this paper we extended the results obtained in [TJ to the case in which the two 
oscillators have different free frequencies, thus when they interact with different parts 
of the reservoirs spectra. We considered the particular case of an Ohmic distribution in 
the high temperature limit for different initial TWB states and different values of the 
free oscillator frequencies. Quantitative differences in the entanglement dynamics can 
be observed in particular when both resonance parameters are in the regime x « 1, 
where the evolution changes also for a small values of Ax. On the contrary the dynamics 
is not strongly affected in the opposite regime even for large detuning. Moreover when 
the two frequencies lay in opposite parameters regime, the dynamics is similar to the 
case x » 1, characterized by lack of oscillations and fast ESD. 
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